Allosteric communication between ACE2 active site and binding interface with SARS-CoV-2.

SARS-CoV-2, the virus causing COVID-19, initiates cell invasion by deploying a receptor binding domain (RBD) to recognize the host transmembrane peptidase angiotensin-converting enzyme 2 (ACE2). Numerous experimental and theoretical studies have adopted high-throughput and structure-guided approaches to (i) understand how the RBD recognizes ACE2, (ii) rationalize, and (iii) predict the effect of viral mutations on the binding affinity. Here, we investigate the allosteric signal triggered by the dissociation of the ACE2-RBD complex. To this end, we construct an Elastic Network Model (ENM), and we use the Structural Perturbation Method (SPM). Our key result is that complex dissociation opens the ACE2 substrate-binding cleft located away from the interface and that fluctuations of the ACE2 binding cleft are facilitated by RBD binding. These and other observations provide a structural and dynamical basis for the influence of SARS-CoV-2 on ACE2 enzymatic activity. In addition, we identify a conserved glycine (G502 in SARS-CoV-2) as a key participant in complex disassembly.

Entropic contribution of ACE2 glycans to RBD binding.

The spike protein of the SARS-CoV-2 virus (the causative agent of COVID-19) recognizes the host cell by binding to the peptidase domain (PD) of the extracellular receptor angiotensin-converting enzyme 2 (ACE2). A variety of carbohydrates could be attached to the six asparagines in the PD, resulting in a heterogeneous population of ACE2 glycoforms. Experiments have shown that the binding affinity of glycosylated and deglycosylated ACE2 to the virus is virtually identical. In most cases, the reduction in glycan size correlates with stronger binding, which suggests that volume exclusion, and hence entropic forces, determine the binding affinity. Here, we quantitatively test the entropy-based hypothesis by developing a lattice model for the complex between ACE2 and the SARS-CoV-2 spike protein receptor-binding domain (RBD). Glycans are treated as branched polymers with only volume exclusion, which we justify using all-atom molecular dynamics simulations in explicit water. We show that the experimentally measured changes in the ACE2-RBD dissociation constants for a variety of engineered ACE2 glycoforms are in reasonable agreement with our theory, thus supporting our hypothesis. However, a quantitative recovery of all the experimental data could require weak attractive interactions.